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Abstract 

The Berezin-Lieb inequalities provide upper and lower bounds for 
a partition function based on phase space integrals that involve the 
Glauber-Sudarshan and Husimi representations, respectively. Gener- 
alizations of these representations have recently been introduced by 
the present authors, and in this article, we extend the use of these 
new representations to develop numerous analogs of the Berezin-Lieb 
inequalities that may offer improved bounds. Several examples illus- 
trate the use of the new inequalities. Although motivated by prob- 
lems in quantum mechanics, these results may also find applications 
in time-frequency analysis, a valuable cross fertilization that has been 
profitably used at various times in the past. 
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1 Introduction 



The Berezin-Lieb inequalities offer upper and lower bounds for partition func- 
tions of elementary quantum systems. In particular, for a system composed 
of a single canonical degree of freedom, let P and Q denote canonical Heisen- 
berg variables, fulfilling the commutation relation [Q, P] = il, in units where 
h = 1. Let |0) denote the normalized ground state of an elementary oscilla- 
tor for which (Q + iP)\0) = 0. Canonical coherent states for this system are 
taken to be states of the form (see, e.g., Refs. [1]) 

\p,q) = U\p,q}\0) , U\p,q] = (1) 

for all (p,q) € K 2 , where U\p,q] denotes the unitary Weyl operator. Let 
H = H(P,Q) denote the Hamiltonian for the system in question. The 
corresponding classical Hamiltonian is denoted by H c i(p,q). We introduce 
two well-known symbols associated with H, namely, the Husimi [2] symbol 
H H (p,q) defined by 

H H (p, q) = (P, q\n(P, Q) \p, q) = {0\H(P + p,Q + q)\0) ; (2) 

and the Glauber-Sudarshan [3] symbol H G _ s (p,q) implicitly defined by the 
operator representation 

H(P, Q) = J H G _ s (p, q) \p, q) (p, q\ dpdq/2n . (3) 

It follows from Eq. (2) that these two symbols are related by the integral 
equation 

H H (p',q') = J \{p',q'\p,q)\ 2 HG-s(P,q)dpdq/27r 

= J e -[(p'-p) 2 +(<7'-<?) 2 ]/2 H G _ s (p, q) dpdq/27r . (4) 
Armed with these definitions, the Berezin-Lieb inequalities [4, 5] read 
J e -PH H (p, q) dpdq/27l < Tr[e-^ p ^} < J e -^ G -s(j>, q ) dpdq / 2ir . (5) 

In what follows we will implicitly rederive this inequality as a special example 
of our generalizations. 
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The purpose of the present paper is to extend such inequalities by offering 
infinitely many additional symbol pairs that can stand in place of the Husimi 
and Glauber-Sudarshan symbols in Eq. (5), thereby generalizing the original 
Berezin-Lieb inequalities. 



2 Multiple phase-space symbols 

In a recent paper [6], the authors have introduced a wide class of phase-space 
symbols that are analogues of the Husimi and Glauber-Sudarshan dual pair. 
Let us first recall the principal elements of that study specialized to the 
discussion at hand. 

We first introduce a nonnegative, trace-class operator a = > which 
we normalize so that Tr(<r) = 1. Such operators have the generic form given 
by 

oo 

(7 = ^d|6,>(6,| , (6) 
i=i 

where denotes a complete orthonormal sets of vectors, and the co- 

efficients {q}^ 1 satisfy the conditions q > and S^q = 1. In short, a 
enjoys all the properties to be a density matrix. 

We shall make use of the function Tr(U[k,x]a) defined for all (k,x) in 
phase space, and we restrict a so that the expression 

Tr(U[k,x]a) ^ (7) 

for all (k,x) G M 2 . 

We next recall the Weyl representation of operators given by 

A = J A(k,x)U[k,x]dkdx/27r , (8) 

where 

A(k,x) = Tr(U[k,x]*A) . (9) 
Given two such operators A and B, it follows that 

Tr(A*B) = J A( y k,x)*B( y k,x)dkdx/2n . (10) 
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In terms of the double Fourier transformation, given by 

A(p, q) = J e i(qk ~ px) A(k, x) dkdx/27c , (11) 

and likewise for B(p, q), it also follows that 

Tr(Atfi) = J A(p, q yB(p, q )dpdq/2n. (12) 

We next modify the symmetric expression for Tr(A^ B) given by Eq. (10) 
so that 

Tr(A^B) = J { T ^ k X) * ]a) } {Tr(U[k, x] a)B(k, x)} dkdx/2n 

= J { Ti ^ k ^ a) V mU[k, x] a)B(k, x) } dkdx/2n 
= J A_ a (k,x)* B a (k,x) dkdx/2-K 

= J A_ a (p, q)* B a (p, q) dpdq/27T . (13) 

In the final line we have introduced the Fourier transform of the symbols in 
the line above. We next show that there are alternative expressions involving 
the symbols A^ a {p, q) and -B CT (j9, q) directly in their own space of definition 
rather than implicitly through a Fourier transformation. 

We begin first with the symbol B a (p, q). In particular, we note that 

B a (p,q) = J e i(fc ^ p) Tr(f/[A;,:r]a)5(A;,x) dkdx/2-K 

= J Tr(U[p, q] f U[k,x]U[p, q]a)Tr(U[k,x] f B)dkdx/2n 

= J Tr(U[k,x]U\p, q]aU[p, q] r ) Tr(U[k, x]^ B) dkdx/2n 

= TT(U[p,q]aU[p,q]^B) , (14) 

where in the second line we have used the Weyl form of the commutation 
relations, and in the last line we have used the Weyl representation Eq. (10), 
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which leads us to the desired expression for B a (p,q). This expression is the 
sought for generalization of the Husimi representation; indeed, if a = |0) (0| 
it follows immediately that 

B a (p,q) = Tr(U\p,q]\0)(0\U\p,q^B) 

= (p,q\B\p,q) = B H (p,q) . (15) 

For general a, to find the expression for A^ a {p, q) we appeal to the relation 

Tr(A*B) = J A_ a (p,q)*B a (p,q)dpdq/27T 

= J A_ a (p,q)*Tr(U\p,q}aU\p,q^B)dpdq/2n , (16) 

an equation, which, thanks to its validity for all suitable operators B, carries 
the important implication that 

A = J A- a (p,q)U\p,q]aU\p,q]Updq/2n . (17) 

Observe that this equation implies a very general operator representation 
as a linear superposition of basic operators given by U[p, q] a U[p, q]^ , for a 
general choice of a. 

Equation (17) for A is the sought for generalization of the Glauber- 
Sudarshan representation; indeed, if a — \ 0) (0|, it follows immediately that 

A = J A_ a (p,q)U\p,q}\0) (0\U[p,q}Updq/2n 
= J A_ a (p,q) \p,q){p,q\dpdq/2n 

= J A G _s(p, q) \p, q) (p, q\ dpdq/2n . (18) 

Once again there is a direct connection between the generalization of the 
Husimi representation, A a (p, q), and the generalization of the Glauber-Sudar- 
shan representation, A_ a (p, q). In particular, it follows that 



A a (r,s) = J A_ a (p, q)Tr(U[r, s]all[r, s]'U[p, q]aU\p, q}*)dpdq/2n 

= J A_ a ( y p,q)Ti(U[r-p,q-s]crU[r-p,q-s] t a)dpdq/27r. (19) 
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This equation is a convolution, which just reflects the multiplicative connec- 
tion between these two symbols in Fourier space. 

3 Derivation of inequalities 

Let {|r)}^2 =1 denote an arbitrary, complete, orthonormal basis. Consider the 
expression [cf., Eq. (6)] 

f(p,q\r) = (r\U[p,q]aU[p,q)}^\r} 

oo 

= J2c l \(r\U[p,q}\b l }\ 2 . (20) 
i=i 

It follows that 

J f(p,q\r)dpdq/27c=l, (21) 

and also that 

oo 

£/(p,?|r) = l. (22) 

r=l 

We can interpret these results in two different ways: On the one hand, 
f(p,q\r) is a probability density on R 2 for each value of r; on the other 
hand, /(p, q\r) forms a discrete probability on {1,2,3,...} for each phase- 
space point (p, q). 

3.1 Jensen's inequality 

The Jensen inequality [7] applies to convex functions <f>(x) — such as e~P x — 
and arbitrary probability distributions on x E K. If ((■)) denotes an average 
over that probability distribution, then the Jensen inequality reads 

< (m) , (23) 

or, in particular, 

e -P(x) < ( e -f)x) _ ( 24 ) 
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This equation will be important in what follows. 

Let T-L denote the Hamiltonian with a discrete spectrum {fj, r }'^ =1 and an 
associated set of eigenvectors { | r) such that 

H\r)=n r \r). (25) 

It also follows that 

oo 

n = J2^\r}(r\. (26) 

r=l 

Following Lieb [5], we first observe that 

(r|e _/W |r) = exp[-fi(r\H\r)] 

= exp[-/3fH_ a (p,q) f(p, q\r) dpdq /2tt] 

< J e'^ H -^ f(p, q\ r) dpdq/2ir . (27) 

Summing on r leads to 

Tr^"^) < J e ^ H -^ q) dpdq/2n . (28) 

Second, we learn that 

exp[-/3H a (p,q)} = exp[-/3£ r /! r /(p,g|r)] 

oo 

< J2 e M-PVr]f(p,q\r). (29) 



r=l 

p2 



Integrating over R leads to 

e -PH a (p, q ) dpdq / 27v < Tr(e^) . (30) 



Above we have two separate inequalities, one an upper bound, the other 
a lower bound. These bounds apply for any choice of a that fits our re- 
quirements, and so we can decouple the choice of a and assert that a can 
be chosen independently in the two cases. In summary, therefore, we have 
established the inequalities 

J e -PH a ,{ P , q ) dpdq / 27r < Tv(e-^ H ) < J e -P«-^i) dpdq/2n , (31) 



where a' and a may be chosen independently of each other. This possibility 
permits optimizing both bounds by taking the supremum over the lower 
bound and taking the infimum over the upper bound. The bounds as given 
by Eq. (31) now lead to upper and lower bounds, respectively, of the free 
energy F(/3) = — In Z(j3)/j3, where Z(/3) denotes the partition function, as 
well as bounds on the ground-state energy Eq since Eq = lim^oo F(j3). 

4 Symbols for the lower bound 

We focus on the symbol 



where the notation (q) a is the symbol H a (p, q) when the operator % is simply 
Q. Below we list a table of symbols needed for our present purposes: 



(q)<t = Q + Q , 
(p)a =P + P , 
(q\ = q* + 2qQ + Q^, 

ip\ =P 2 + 2pP + P 2 , 

(qp)a = qp + qP + pQ + QP , 

(pq)a = pq + pQ + qP + PQ , 

(q% = q 4 + 4g 3 Q + Qq 2 W + W + Q 1 , 

(p% =P 4 + 4p 3 P + 6p 2 P2 + 4 p p3 + pi ; 

(q 2 p 2 )a = qY + 2pq 2 P + 2qp 2 Q + q 2 P^ + pW + MpQP 

+2qQP 2 + IpWP + WP 2 , 
(P 2 q 2 )a = P 2 q 2 + 2pq 2 P + 2qp 2 Q + q 2 P* + pW + ^pqPQ 

+2qP^Q + 2pPQ^+P^ . (35) 




(32) 



(33) 



In that case we find, e.g., that 



(q) a = TT((Q + q)a)=q + Q , 



(34) 
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Note that on the left-hand side the order matters, i.e., (qp) a ^ (poja, etc. 
We also notice that for the quadratic symbols 

(<? 2 V = (g + Q) 2 + A(Q), 

(p 2 ) CT = (p + P) 2 + A(Ph 

\[{qp)* + ipq)a] = (q + Q)(p + P) + A(Q, P) , (36) 

mterms of the variances A(0) = O^-Q 2 and A(O u 2 ) = {C\0 2 + 2 0{)/2- 
0\ 2 . For a conventional minimal uncertainty state, e.g., A(Q)A(P) = 1/4 
and A(Q,P) = 0. 

We also introduce a special-case table based on a symmetry we shall 
impose on a, and to be made use of below, namely, that all odd-order averages 
vanish, i.e., Q = Q 3 = P = P 3 = Q 2 P = 0, etc. This special-case table reads 

(q)a = q , 

(q\ = q 2 + W, 

(p\ = p 2 + p\ 

(qp)a = qp + QP , 
(pq)a = pq + PQ , 
(q% = q 4 + Qq 2 W + W, 
(p 4 )a = p A + Qp 2 P* + P I , 

{q 2 p 2 )a = qV + q 2 P* + P 2 W + MpQp + Wp 2 , 

{P 2 q 2 )a = pV + q 2 P^ + pW + 4pqPQ + PKp . (37) 



5 Symbols for the upper bound 

The construction of the upper limit is somewhat more involved than that for 
the lower limit. We start with Eq. (17), which is 

A = J A_ a (p,q)U[p,q}aU\p,q}Updq/2K . (38) 

For reasons of clarity we limit ourselves to a a number-operator diagonal form 
for a, i.e. a = S^L c n |n)(n|, c n > 0, and S^L c n = 1, where N\n) = n\n) 
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for the number eigenstates {|n)}. We learn that in general 

A=^ =0 c n J A_ a (p,q)U[p,q]\n)(n\U[p,q}Updq/2n 

= ^ =0 c n J A^ a (p,q)\p,q;n)(p,q;n\dpdq/2ir , (39) 

in terms of the so called semi-coherent states or displaced coherent states 
\p, q;n) = U[p,q]\n) (see, e.g., [8]). To see what this means, let us take a 
simple example with A = P 2 + Q 2 . Since an operator is determined by its 
expectation value in canonical coherent states, it is sufficient to consider the 
Husimi symbol An(p,q) as given by Eq. (2), i.e., 

(r, s; 0| (P 2 + Q 2 ) \r, s; 0) = (0| [(P + r) 2 + (Q + s) 2 ] |0> = (r 2 + s 2 ) + 1 

= S^ =0 c n J [(p + r) 2 + (q + s) 2 + k 2 ] p n (p, q) dpdq/2^ , (40) 

where 

Pn(p,q) = |(0|p,g;n)| 2 = e -^ + ^ 2 {p 2 + g 2 )«/2"n! . (41) 
Here, we have made use of the Ansatz 

(P 2 + Q 2 )_ (T (p,q)=p 2 + q 2 + k 2 , (42) 
where k 2 is a constant to be determined, and we immediately learn that 

k 2 = -1 - 2E~ =0 c n n = -1 - 2n , (43) 



where we have defined mean-values f(n) = S^ =0 c n /(n). It now, e.g., follows 
that the right-hand side of Eq. (17), with the upper symbol as given by 
Eqs. (42) and (43), has \n) as an eigenstate with eigenvalue 2n + 1. It is not 
entirely trivial to verify this explicitly, but it follows using the properties of 
displaced coherent states as well as properties of the conventional associated 
Laguerre polynomials L™\ 



k=0 
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In like fashion, it follows for A = (P 2 +Q 2 ) 2 and the corresponding Husimi 
symbol that 

(r, s; 0| (P 2 + Q 2 ) 2 |r, s; 0) = (0| [(P + r) 2 + (Q + s) 2 ] 2 |0) = (r 2 + s 2 + l) 2 
ee ^ =0 c n | [((p + r) 2 + (q + S ) 2 ) 2 + fc 4 ((p + r) 2 + (g + S ) 2 ) + k 6 ] 

xp n (p,q)dpdq/2n , (45) 

expressed in terms of the (assumed) symbol 

((P 2 + Q 2 ) 2 )-^, q) = (p 2 + q 2 ) 2 + Hp 2 + g 2 ) + k 6 . (46) 

One now finds, making use of Eq. (41), that 

A; 4 = 2-8£^ =0 c n (n + l) = -6-8n, (47) 

and 

fc 6 = 1 - 4E~ c n (n + l)(n + 2) - 2£; 4 E~ c n (n + 1) 

= 5 + 16n + 16n 2 - 4^ . (48) 

In a similar manner and for A = Q A , we can write 

(Q 4 )_ CT (p,g) = g 4 + a 2 g 2 + a 4 , (49) 

where 

a 2 = -3(1 + 2n) , (50) 

and 

a 4 = 3(I + |n + 2n 2 -i^) . (51) 

The expressions above now relate the standard symbols to the generalized 
symbols. Extension of these expressions to other polynomials in P and Q is 
straightforward. 
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6 Examples 



With the special choice for a considered above, i.e., a — ££L c n |n)(n|, c n > 0, 
we will now consider some specific examples in order to 
illustrate the use of the generalized upper and lower symbols. We first remark 
that in the trivial case of an harmonic oscillator with H = (P 2 + Q 2 )/2, such 
that Z(f3) = l/[2 sinh(/3/2)], the lower symbol Eq. (36) and the upper symbol 
Eq. (42), together with the bounds Eq. (31), lead to the expression 

e -/3(A(P)+A(Q))/2 < Z ((3)< e / J ( 1 /2+n)/2 ( 52 ) 

which, obviously, is true. We can optimize this expression in the form 

e-^ 2 //? < Z(f3) < e V 2 /p . (53) 

From the corresponding lower bound we then obtain an upper bound on the 
ground-state energy E < 1/2 since E = — lim^^oo In Z(f3)/(3. In the high- 
temperature limit, i.e., /3 — > 0, the bounds in Eq. (53) exactly reproduce the 
classical Gibbs partition function Z d (f3)/2ii — 1/(3 taking the fundamental 
phase-space volume 2n into account and making use of 

Z d {(3) = j e-^^dpdq, (54) 

with, of course, H d (p, q) = {p 2 + q 2 )/2. 

6.1 A non-linear oscillator 

Here we consider Hamiltonians of the form % = "H(iV), where N is the usual 
number operator. We study this example more for its ease of analysis and 
pedagogical value. We choose as our example H — (N — a)(N — b). Such a 
form of an Hamiltonian has its roots in, e.g., the description of a single- mode 
non-linear Kerr-medium in quantum optics or a single vibrational mode be- 
yond the harmonic approximation. We make the choice a = 1 and 6 = 5. 
We observe that the partition function Z(j3) = J2™=o ex P[— P( n ~ ^)( n ~ 5)] 
then has the form Z(f3) ~ exp(A/3) for large values of (5. A straightforward 
application of Poisson re-summation techniques also leads to the behavior 
Z(f3) ~ y/nJP/2 for small values of f3, which corresponds to the high- 
temperature limit of the classical partition function Z c i/2ti using Eq. (54) 
with H d = (p 2 + q 2 f/A - 7(p 2 + q 2 )/2 + 33/4. 
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We may then combine these factors for TL = (N — 1)(N — 5) at hand by 
noting that 

(AT - i)(N - 5) = \ (P 2 + Q 2 - l) 2 - 6 \ (P 2 + Q 2 - 1) + 5 

= i(P 4 + Q 4 + P 2 Q 2 + Q 2 P 2 ) - 7±(P 2 + Q 2 ) + 33/4 . (55) 

Consequently, 

Ha(p, q) = \[{p% + (q% + (p 2 q\ + (qVU - 7±[(p 2 ) CT + (g 2 ),] + 33/4 
= |[p 4 + Qp 2 ! 52 + P 4 + g 4 + Gg 2 ^ 2 + Q 4 + g 2 p 2 + q 2 ! 52 + p 2 ^ 2 
+4qpQP + Q^P 2 + p 2 q 2 + q 2 ! 52 + p 2 ^ 2 + ApqPQ + P~ 2 Q 2 } 

-7|[p 2 + i^ + g 2 + g I ]+33/4. (56) 

Since we have restricted our choice of er so that it is only a function of 
N, i.e., a = cr(N), a has now a symmetry that makes Q 2 = P 2 = C 2 , 
p4 = gl = c 4) g2p2 = p2g2 = ^ and importantly that ~QP + P~Q = 0. 
The three constants C 2 , C 4 , C 22 are the only remnant of a in H a (p,q), and 
of necessity, they satisfy C 2 > 1/2, C 4 > Cf, and C 4 > C 22 . With the 
restriction cr = cr(N) we can actually be more precise and write 

C 2 = \+n, C 22 = I(^ + n+i) , C 4 = |(^ + n + i). (57) 

Putting this information together, we find that 

H„(p, 0) = W + q 2 f + K x {p 2 + q 2 ) + K 2 , (58) 

where 

K 1 = l(C 2 -2) = l(n-l), 

K 2 = C 4 + \C 22 - 7C 2 + f = (n - 3) 2 - 4 . (59) 

We note the fact that H a (p, q) is a function only of the combination 
(p 2 + g 2 ) on the basis of our restriction that a = cr(N). It follows, therefore, 
that the lower bound of interest is given by 

J exp{-PH a (p,q)} dpdq/27r = If™ exp{-Pl\s 2 + K lS + K 2 ]}ds , (60) 
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where we have passed to polar coordinates and set s = (p 2 + q 2 ). The upper 
bound integral is a function of j3 as well as the a-parameters, C 2 , C4, and 
C 2 2, i.e., the independent mean- value n and dispersion (n — n) 2 parameters. 

The lower bound of Eq. (31) together with Eq. (60) now leads to the lower 
bound ^/"7r//3/2 < Z(f3) as f3 — > 0. This lower bound again corresponds to 
the high-temperature limit for the classical partition function Z c i(/3)/2 / k. By 
making use of E = — lim^^oo In Z(j3)/j3, Eq. (60) leads to the upper limit 
Eq < —4 using the state a = |3)(3|. We observe that such a state will not 
strictly satisfy the restriction imposed by Eq. (7) since Tr(U[k,x]a) then 
will be zero at isolated points away from the origin k — x — 0. But, in fact, 
the restriction Eq. (7) is then not required if A is a polynomial in P and 
Q since the symbol A(k,x) as defined in Eq. (9) will involve derivatives of 
delta-functions with support at the origin [9, 6]. 

The upper bound of Eq. (31), using Eqs. (42) and (46), now leads to 

Z(P) < \ / °° exp{-/3[|(s 2 + k,s + k 6 )-l(s + k 2 ) + f]}ds , (61) 

where the parameters k 2 , k^ and ke are given by the equations (43), (47) 
and (48), respectively. It is now evident again that Eq. (61) reproduces the 
high-temperature limit of the classical partition function Z c i(f3) ~ a/vt//3/2. 
The upper bound of Eq. (31) gives unfortunately now a rather poor lower 
bound on the ground state energy E > —12 — 9n — n 2 , i.e. E > —12. 



6.2 An anharmonic oscillator 

We next consider the Hamiltonian H = (P 2 + Q 2 )/2 + \Q 4 /2 > 0, A > 0, to 
define the partition function. With the lower and upper symbols as given by 
Eqs. (37), (42), and (49), we now find that 



Z(0) < 1 c ~/3fe+Aa 4 )/2 f e ~(x 2 +\(x*/p+a 2 x 2 ))/2 dx 

' ~ pV2^ J 



(62) 



and 



Z(B)> 1 c -/3(A(P)+A(Q)+AQ4)/ 2 f e -(x*+\{x± / p+^Q 2 ))/2 d% /ggx 

~ PV2^ J 

In the limit of large /3, the lower bound on Z(/3) and the fact that H > then 
lead to < E < (1 + \Q I )/2. With a = |0)(0| one finds the upper bound 
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Eq < (1 + 3A/4)/2 which, e.g., can be compared to the "exact" numerical 
value of 2E = 1.392351641530... for A = 1 [10]. We expect that this upper 
bound could be improved with a different choice of a. 

A consequence of the upper and lower bounds Eqs. (62) and (63) now is 
that for sufficiently small (3 the upper and lower bounds converge to the well 
studied (see, e.g., Refs. [11] ) classical and asymptotic form 

Z{fi) = Jj^J e-^'^^dx ee Z d {P)/2* 

^ 16X K 1/4 ^/l6\) (64) 



2Av / 2¥V ^ 

using Eq. (54) with H c i(p, q) = (p 2 +g 2 )/2 + Ag 4 /2. The expression in Eq. (64) 
involves all the energy states of the anharmonic oscillator in a highly non- 
trivial manner. In our case we are specifically interested in the limit /3 — > 0, 
i.e., Z(/3) ~ r(l/4)(2/3/A) 1 / 4 /2/3 v / 2T. 



7 Comments 

For clarity, we have mainly focused on matrices a = cr(N) which meant that 
a = S^l c n |n) (n\. More general matrices of course would involve expansions 
of the form 

(T = ^n,n'=O c n,n'\n){n'\ (65) 

expressed in terms of a general matrix {c„ jn /} that still ensures that a has all 
the properties of a partition function. The use of such more general choices 
for a will inevitably lead to expressions involving the matrix elements [8] 

(n\U\p,q]\n') 

= \J^M-\{P 2 + <f)]{<l + ^P) n - n ^-'(|(P 2 + g 2 )) , (66) 

for n > n' expressed in terms of the associated Laguerre polynomials Eq. (44); 
instead, when n < n', use (n| U\p, q] \n') = (n'\U[—p,—q]\n)*. The simple 
example where H = P 2 + u 2 Q 2 , u ^ 1, shows that the optimal choice of 
a is not always given by |0)(0|, where (Q + iP) |0) = 0, but in the present 
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case by a = |0; oo)(0; u\, where (uQ + iP)\0;u) = 0. This remark serves 
to confirm that the generalized representations have the possibility to make 
better bounds. It may be true that choices for a of the form (analogues 
of pure states) may be optimal, and that perhaps choosing as the ground 
state of the Hamiltonian under examination may lead to optimal bounds. 
Those are interesting questions for the future. 

8 Conclusion 

We have developed new, classical, phase space bounds to deal with special- 
ized (i.e., the partition function) questions that arise in quantum mechanics, 
and which, by their very nature, are technically easier to deal with than in 
their original form. It is quite likely that the generalized phase-space sym- 
bols we have introduced may have additional applications both in quantum 
mechanics and in time-frequency analysis. 
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